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Abstract. 

We report results from a high-resolution optical spectroscopic survey aimed to search for nearby young associations 
and young stars among optical counterparts of ROSAT All-Sky Survey X-ray sources in the Southern Hemisphere. 
We selected 1953 late-type {B — V > 0.6), potentially young, optical counterparts out of a total of 9574 IRXS 
sources for follow-up observations. At least one high-resolution spectrum was obtained for each of 1511 targets. 
This paper is the first in a series presenting the results of the SACY survey. Here we describe our sample and our 
observations. We describe a convergence method in the (UVW) velocity space to find associations. As an example, 
we discuss the validity of this method in the framework of the /3 Pic Association. 
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1. Introduction 

Post-T Tauri Stars (pTTS), and its prototypical exam- 
ple FK Ser llHerbid Il97.'^. were intr oduced as a new, 
systematic, category bv lHerbi'3 (^7^ as T Tauri stars 
(TTS) that follow classical and weak line TTS (cTTS and 
wTTS) in an evolutionary sequence. The pTTS were ex- 
pected to outnumber the cTTS and wTTS populations by 
far, but their discovery and identification remained diffi- 
cult f or a long time. Th e relation of pTTS to "isolated" 
TTS l|Quast et al.lll987(l (young, low-mass pre-main se- 
quence stars found spatially far away from any apparent 
dark or parental molecular cloud) remained an open prob- 
lem. Interestingly, Herbig's catalog of pTTS, and the sug- 
gested list of isolated TTS produced independently by the 
Brazilian group are actually quite similar. An intriguing 
case is T W Hya, a high Galactic lati tude object, confirmed 
as TTS (|Rucinski fc Krautteilll983|) . with a distance of at 



* Based on observations made under the ON-ESO agreement 
for the joint operation of the 1.52 m ESO telescope, on two 
2.2 m ESO telescope runs, and at the Observatorio do Pico dos 
Dias MCT/Laboratorio Nacional de Astroffsica (LNA/MCT), 
Brazil 

** Tables 3, 4 and 5 are only available in electronic form at the 
CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or 
via,http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A-|-A/ 



least 13° from the nearest dark clouds. Its origin remained 
enigmatic for a long time. 

A systematic search for more isolated TTS was pur- 
sued with the optical spectroscopic Pico dos Dias Survey 
(PDS) among optical counterparts of the IRAS Point 
Source Catalog JCregorio-Hetem et al ] ll992t iTorres et alJ 
1199,4 lTorresl1l999l) . One of the first results of the 
PDS was the discovery of four additional TTS around 
TW Hya l)de la Reza et al.l Il989l: iGreeorio-Hetem et alJ 
1 1 992fi . They concluded that this group was likely a very 
young association relatively close to the sun. 

The putative TW Hya association was actually a very 
odd resTilt. In fact, in a recent review. ^Zucke rman & Sond 
(2004) reminded us that until the late 1990 the Hyades 
and the UMa clusters (with an age about 600 Myr and 
300 Myr, respectively) were the only coeval, co-moving 
concentrations of stars known within 60 pc from Earth. 
Thus, the closest place to find very young stars (typically 
about 10 Myr) was in the nearby classical star forming 
regions, such as Sco-Cen, Taurus, p Oph (at ~150 pc) 
and Orio n somewhat far t her at ^450 pc. 

Later, iKastner et al ] ^ll997^ confirmed that TW Hya 
forms a physical association of about 20 Myr and at a dis- 
tance of 40 to 60 pc from Earth - the TW Hya Association 
(TWA), based on the similarity of the X-ray fluxes, radial 
velocities, astrometry (Hipparcos) and their spectroscopic 
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characteristics. Even additional members were found later 



llWebb et al.l llQQfi 


ISterzik et al.l 119991: 


l2001aHSons- et all 





Only very few good candidates for isolated TTS were 
found with the PDS. But while the IR-excess selection 
criteria effectively finds young stellar objects embedded in 
their placental material or with circumstellar disks, it fails 
to signal older objects whose disks have been dissipated. 
Hence most stars with ages between 10 — 70 Myr (i.e. 
wTTS and pTTS) escaped the discovery by this method. 

Due to th e enhanced X-ray activity of young stars 
llWalteilll986t) more efficient selection criteria for post and 
isolated TTS candidates were developed. The high sen- 
sitivi ty and full sky coverage of the ROSAT all-sky sur- 
vey l)TruemDeilll982|) revealed thousands of new X-rays 
sources pro jected in the directi on of nearby star form- 
ing regions (|Guillout et al.l[l998|) . Ground-based spectro- 
scopic studies showed that a large fraction of these X-rays 
sources were indee d wTTS together with older pTTS and 
ZAMS stars, fe.g. lAlcala et al. 2 0Q0j . Surprisingly, many 
of the newly found weak-line TTS were not obviously con- 
nected to any molecular cloud regi on which again raised 
many questions about their origin l)Sterzik ^t al. 199,^/1 . 

Young, nearby associations, and in particular TWA, 
are of great importance to understand the local star form- 
ing history, and they reveal samples of stars that al- 
low to study the transition phase of disk dispersal and 
planet formation in great detail. Therefore several groups 
have started to look for similar nearby young associa- 
tions hidden among the ROSAT X-ray sources. Shortly 
after the discovery of TWA eight mo re young nea. r by as- 
sociations have been announ ced fsee iTorres et el\ l(2003j) 
and the review by Zuckcr man fc Sond 2004l. namely, 
77 Cha (d ~ 97 pc a nd age ~ 9 Myr; Mamaiek et al.l 
Il999t iLvo et all l2004l): e Cha id ~ 115 pc and age 
~ 5 Mvr: lMamaiek et aLlExxil iFeigelson et al.ll2n0.ll) : the 
Horologium Associ ation (HorA) (d ~ 60 pc; age ~ 30 Myr; 
iTorres et al.l l2000l) and the adjacent and very similar 
Tucana Associati o n (Tuc A) simultaneously discovered by 
IZuckerman et al. I ll2001bl) : th e Pic Association (BPA 
id ^ 35 pc; age ^ 15 Myr; iBarrado v Navascuesllm 
IZuckerman et al.ll2001a^; the AB Dor Association {d ^ 
20 pc ; age ~ 50 Mvr: IZuckerman et alJEooi iTorres et alJ 
I2OO3I - calling it as AnA); and finally two less well de- 
fined asso ciations, the Argus Association {d ^ 100 pc; age 
~ 30 Mvr: {Torres et a L 20031 and the Octans A ssociation 
(d - 110 pc; age ~ 30 Myr; Torres et alJl2003l) . 

The similarity between the adjacent TucA and HorA 
motivated us to start SACY (Search for Associations 
Containing Young stars) ^ survey in 2000. Some early ob- 
sevations, obtained in 1999 for other purposes, were also 
included in the SACY, whose initial goal was to exam- 
ine the physical relation between these two associations. 
Later, SACY was enlarged in order to look for other, hith- 



^ Sacy is a mythological Brazilian being, a little black boy 
with only one leg, that frightens night travelers. 



erto undiscovered, nearby young associations among the 
optical counterparts of the ROSAT X-ray sources. 

Preliminary results of the SACY h ave been presented 
in different conference contributions ijTorres et alJl2003. 
and references therein), where we have suggested eleven 
potential nearby associations using the preliminary data 
(including components of the Sco-Cen Association). One 
of those first results was the existence of two some- 
what superposed associations which we called GAYAl and 
G AYA2-^ . GAYAl is so si milar to the Tuca na/Horologium 
Association proposed bv IZuckerman fc So ng ( that 
we can adopt this designation, although some of their 
proposed members seem to match better with GAYA2 
[d ~ 80 pc; age ~ 20 Myr), as we will present in a forth- 
coming paper. In fact, the present paper is the first one 
of a series that aims to present the results of the SACY. 
Here, we present the survey, the observed sample and the 
observing strategy (Sec. 2). In the associated catalog, we 
give the spectroscopic properties (radial velocity, Fsin(i), 
spectral type. Lithium 6708 A and Ha equivalent widths) 
as well as the UBV{RI)c photometry of the observed 
stars. Our kinematical appr oach is presented and ex em- 
plified with the known BPA jZuckerman et aLll2001a[) for 
which we propose new members. 

2. The sample and observations 

The ROSAT All-Sky Bright Source Catalogue (IRXS) 
ijVoges et al.l [l999) was used as a starting point to 
build an initial sample containing young star candidates. 
Therefore, IRXS was cross-correlated with the Hipparcos 
and Tycho-2 Catalogue s. We followed an approach simi- 
lar to that described bv lCuillout et all §999). Within an 
error radius of 2.6 times the ROSAT positional error one 
(or more) possible optical counterparts were selected for 
each IRXS source. Only stars having B ~ V > 0.6 (from 
Hipparcos or Tycho-2) were considered. This cut-off is jus- 
tified by the lithium depletion behavior. For young stars 
with spectral types later than mid-G (i.e., B — V ^ 0.7) 
the strength of the Li 6708 A line can be used as an youth 
indicator (Martin 1997). Our cut-off (0.6) corresponds ap- 
proximately to a GO dwarf, in order not to loose interesting 
objects close to this boundary. We exclude hotter stars, be- 
cause they lack of a similar youth indicator and due to the 
difficulty to obtain precise radial velocities (Vr) necessary 
for the convergence method later employed in SACY. 

We also excluded all Hipparcos stars having 
A'lv < 2.0, because for our color limit young stars older 
than 5 Myr must be less luminous than IOLq. There is 
some subjectivity in our final list, as the positions, veloci- 
ties, parallaxes, colors, and magnitudes are all affected by 
different errors, and we had to decide on a case by case 
basis about inclusion in our final list. For example, some 
late F stars were included and some early G may have 
been excluded, but these cases are not very important as 
we use a conservative color limit. 



GAYA stands for Great Austral Young Association. 
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Fig. 1. Celestial polar projection of the stars observed in SACY. Stars classified as young according to their Li I 
equivalent width (see section 4.2) are plotted as filled circles. More evolved stars are indicated as open circles. The 
transverse curve represents the Galactic plane.) 



Table 1. Completeness of SACY. The table gives the 
northern boundary declination as a function of right as- 
cension up to which the SACY is complete. 



RA interval 


Declination 


(h) 


n 


01-05 


-20 


05-09 


-25 


09-16 


-20 


16 - 18 


-15 


18 - 22 


-10 


22 - 01 


+10 



From an initial list containing 9574 IRXS sources in 
the Southern Hemisphere, we found 1953 of them having 
counterparts in the Tycho-2 and Hipparcos Catalogues us- 
ing the above criteria. During the last seven years, 1511 



stars were observed at least once. Data from the litera- 
ture for 115 objects were used to complement the observed 
sample. (Other 153 objects were observed that do not ful- 
fill the SACY criteria.) 



In Figure^we show the observed sample and in Tabled 
we give the completeness of SACY as a function of right as- 
cension. Only 3 selected sources were not observed within 
this area and there are also 29 stars observed outside this 
area, mainly near the limits defined in the table. Some 
of the holes seen in Figure ^ may be due to incomplete- 
ness of IRXS. The concentration of young stars between 
RA ^ 12 h to ~ 17 h, corresponding to the Sco-Cen 
Association, is clearly visible in Figure Q 
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Table 2. The rejected sample. Column 1: ROSAT 
Identification; Column 2: Identification from Hipparcos or 
Tycho-2 Catalogues; Column 3: Spectral type or other rea- 
son for rejection; Column 4: Other identification. 



IRXS 


Idcnt 


Type 


Other 


(1) 


(2) 


(3) 


(4) 


J005307. 5-743903 


Hip 4157 


RS 


CF Tuc 


J010945. 2-201300 


Hip 5452 


WUMa 


CT Get 


J011143. 2-255720 


T6425-1882 


KOIII+ 


BZ Scl 


J014852. 7-205340 


Hip 8447 


WUMa 


TW Cot 


J015951. 5-380725 


Hip 9336 


K2III 


HD 12320 


J023401. 2-653634 


Hip 11934 


WUMa 


WY Hor 


J024208. 5-463121 


Hip 12611 


KlIII 


HD 17006 


J024326. 2-375540 


Hip 12716 


RS 


UX For 


J035941. 3-584028 


Hip 18659 


KOIII 


UY Ret 


J040729. 6-523413 


Hip 19248 


RS 


AG Dor 


J041425. 3-204921 


T5889-1407 


KlIII SB2 


HD 26917 


J044145. 7-471540 


Hip 21843 


G8III 


HD 30042 


J050135. 3-444955 


Hip 23381 


KOIII 


HD 32517 


J051833. 2-681328 


T9162-0347 


WUMa 


RW Dor 


J082401. 6-694149 


T9198-0008 


A2 


HD311345 


J085943. 2-274850 


Hip 44164 


RS 


TY Pyx 


J091756. 4-372350 


Hip 45623 


G3III 


HD 80332 


J100352. 1-561134 


T8603-1735 


KlIII 


HD 87525 


J103733. 8-564757 


T8609-1385 


XB 


4U1036-56 


J110952. 1-762909 


T9410-2627 


A7V 


CD-75 522 


J111510. 2-611540 


Hip 54948 


CI 


NGC3603 


J112545. 4-763035 


T9411-1577 


KOIII 


HD 99558 


J112557. 8-401547 


T7743-1091 


G6II 


HD 99409 


J113616. 1-380210 


T7740-0698 


RS 


V858 Con 


J113922. 1-392315 


Hip 56851 


RS 


V829 Con 


J114929. 5-504849 


T8228-0126 


KlIII 


HD102726 


J125850. 0-365827 


Hip 63347 


WUMa 


V839 Cen 


J132711. 7-245126 


T6717-0230 


Gal Sey 




J135155. 7-363726 


Hip 67682 


WUMa 


V757 Con 


J140710. 6-510131 


T8276-0897 


Gal 




J160546. 1-395032 


T7855-0905 


D2"(A5V-I-?) 


A:HD144118 


J162804. 7-491150 


T8320-2046 


XB 


4U1624-49 


J163528. 5-480556 


T8329-1465 


B8 


HD330993 


J172446. 8-341224 


T7383-0007 


CI 


Pr 24 


J172509. 3-341117 


T7383-0162 


CI 


Pr 24 


J173728. 0-290759 


T6839-0472 


Gal Sey 




J175620. 0-721850 


T9297-0956 


F2IV DO. 8" 


HD162128 


J180108. 7-250444 


T6846-0633 


XB 


4U1758-25 


J180847. 7-260037 


T6847-3167 


B9 


HD315209 


J181806. 3-121439 


Hip 89681 


CI 


NGC 6604 


J181842. 3-134706 


T5689-1381 


CI 


NGC 6611 


J182030. 8-161022 


T6265-1255 


CI 


NGC 6618 


J182615. 1-145034 


T5702-1197 


XB 07V 


V4138 Sgr 


J192240. 0-203840 


Hip 95244 


K2III 


V4138 Sgr 


J202052. 5-300213 


T7442-0726 


Gal CI 


AGO 3674 


J213949. 1-160018 


Hip 106961 


RS 


AD Cap 


J214906. 4-304157 


T7488-0390 


Gal CI 


AGO 3814 


J220036. 3-024433 


Hip 108644 


G5IIH-sdOB 


FF Aqr 


J231324. 1-1-024028 


Hip 114639 


KlIII 


SZ Psc 


J232748. 5+045126 


Hip 115819 


WUMa 


VZ Psc 


J234351. 0-151655 


Hip 117054 


M4III 


R Aqr 



2.1. The rejected sample 

Whenever available, SIMBAD information was used to ex- 
clude stellar counterparts of the 1953 southern sources 
being well known RS CVn, W UMa, giant stars or any 
other contamination showing ah initio clear signs that 
they are not young late type stars. Within the area defined 
in Tabled we rejected 51 sources and they are presented 
in Table El 

In such way, we have covered ~ 80% of the 1953 
sources. 



2.2. Spectroscopic observations 

Most of the spectroscopic observations 70%) were 
performed with the high-resolution (i? ~ 50000) 



FEROS echelle spectrograph jKaufer et al.lll999|) at the 
1.5 m/ESO telescope at La Silla (Chile) between January 
1999 and September 2002 (ON-ESO agreement and 
ESQ program identification 67.C-0123). In October 2002, 
FEROS was moved to the 2.2 m telescope. After that, 
two more runs (ESO program identifications 072.C-0393, 
077.C-0138) were carried out aimed mainly to collect ad- 
ditional radial velocity points to resolve the orbits of the 
previously identified spectroscopic binaries. We noticed no 
relevant changes in the zero-point and/or in the spectro- 
graph resolution (see below). The FEROS observations 
were taken in the OS (Object-Sky) mode. The reduction 
of the spectra was performed by the FEROS pipeline in- 
cluding flat-fielding, background subtraction, removal of 
cosmic rays, wavelength calibration and barycentric veloc- 
ity correction, yielding as the final product a ID re-binned 
spectrum. 

Another set of data 30%) was collected at the 
coude spectrograph attached to 1.60 m telescope at the 
Observatorio do Pico dos Dias (OPD), LNA, Brazil. The 
setup used was centered at 6500 A covering a spectral 
range of 450 A with a resolution of i? ~ 9000. These ob- 
servations were reduced with standard IRAF packages. 

A few observations were collected using the two- 
fiber-fed hi£h;Tesolution_spectrograph CORALIE [R ^ 
47000) l|Queloz et al.l 12000^ attached to the Swiss Euler 
Telescope, also at La Silla, Chile. As for FEROS, all ob- 
servations were taken in the 0BJ2 mode, i.e., one fiber 
centered on the target star and the other fiber illuminated 
by the background sky. Similarly to FEROS, CORALIE 
also has an on-line reduction system which processes the 
raw frames following the same steps as described above 
for FEROS pipehne. 

2.3. Photometric data 

U BV{RI)c photometry f or part of the sample was ob- 
tained using FOTRAP l|.Tablonski et all Il994l) at the 
0.60 m Zeiss telescope of the OPD. When a star was 
not observed photometrically by us we try to obtain 
some useful photometric data in the Hipparcos and Tycho 
Catalogues or in the available literature in the SIMBAD. 
For multiple stars, magnitudes and colors were corrected 
in order to account for the presence of the companion. We 
will use in the paper preferentially the {V — I)c color. If 
no observed value is available we use one deduced from 
the {B — V)t when F < 10 or otherwise from the spectral 
type. 

3. The SACY Catalog 

The SACY Catalog contains the objects studied in our 
survey. The tables are only available in electronic form. 
The Table 3 (The SACY sample) contains the data for 
1626 stars, 1511 observed by us and 115 objects with in- 
formation from the literature. In the Table 4 (Other ob- 
served stars) we hst 153 observed stars that do not 
fulfill the SACY criteria. The structure of both tables is 
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identical, and they contain astrometrical data, radial and 
rotational velocities, photometric and spectroscopic data. 
Bibliographic sources used in the catalog are given in the 
Table 5 (The bibliographic sources). 

4. Analysis 

The primary goal of the SACY is to derive space motions 
for all young IRXS sources in the Southern Hemisphere. 
Thus, our main effort was to obtain reliable radial veloc- 
ities from the spectroscopic observations. Combined with 
Hipparcos parallaxes and proper motions, they allow us 
to compute the (UVW) components of the Galactic space 
velocity vector which, in turn, is the key ingredient to find 
young associations (see below). The spectroscopic obser- 
vations also enable us to obtain the spectral types and to 
measure both the Li 6708 A and the equivalent widths, 
which are used as youth indicators. 

4.1. Radial and rotational velocities 

Radial velocities were derived by cross-correlating our 
FEROS and coude spectra with the spectrum of the ra- 
dial velocity standard Tau Cet (HD 10700) observed at 
the same night and used as template. During the peri- 
ods where Tau Cet was not visible, 3 other radial velocity 
standards from the CORALIE Extra-Solar planet Survey 
were kindly provided by the Swiss Team and used as ref- 
erence. In very few nights no template could be observed. 
In these unfortunate cases, the template observed in an 
adjacent night was used. As we have not used very red 
templates, our radial velocities for red dwarfs may be less 
precise. 

In order to spot eventual errors, radial velocities were 
computed by two of us independently using the task fx- 
cor fro m IRAF and the method described in Melo et aD 
l|200l|^ . Thi s latter makes use o f KO CORAVEL-type bi- 
nary mask l|Baranne et aljfl979(l as cross-correlation tem- 
plate. In both methods, the resulting cross-correlation 
function (CCF) can be approximated by a Gaussian func- 
tion whose center readily gives the radial velocity and the 
width is related to the broadening mechanism such as tur- 
bulent motions, gravity pressure and rotation. 

The photon noise errors on th e radial velocity m ea- 
surements computed as described in lMelo et al.l l)200l|) are 
very low thanks to our good S/N which is always better 
than 30-50. Thus the uncertainties on the radial velocities 
are dominated by the spectrograph drift during the night 
due to changes in the air index and atmospheric pressure. 
These shifts are of the order of a few hundreds m s~^ per 
night. A more realistic idea of our overall precision can be 
obtained by measuring the dispersion of the radial velocity 
standard Tau Cet measured almost every night during our 
observing runs. The mean radial velocity and its r.m.s are 
respectively -16.316 km s"^ and 0.320 km s^^ spanning 
over 1050 days of interval. 

In Figure El (top) we com pare our radial velocity 
measurements with those from iNordstrom et all l|^004) . 
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Fig. 2. Comparison of our derived radial velocities and 
ysin(z) with those derived bv Nordstrom et al.. (,200^ . 
Top. Radial velocities obtained with FEROS are shown 
as filled circles whereas those obtained with the coude are 
seen as filled squares. The continuous line is the one-to- 
one relation and the two dashed-dotted lines delimit a 
±5 km s~^ region around the continuous line. Bottom. 
Again, the one-to-one relation is shown as a continuous 
line whereas the dashed-dotted lines mark the V sin(i) er- 
rors of 2 km for T^sin z < 20 km s^^ and 10% of the 
observed V sa\{i) for V sni i > 20 km s^^. Multiple stars 
are shown as encircled symbols. 

Radial velocities obtained with FEROS arc shown as 
filled circles whereas those obtained with the coude are 
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seen as filled squares. The continuous line is the one-to- 
one relation and the two dashed-dotted lines delimit a 
±5 km region around the continuous line. We have 
used the mean radial velocity error (e(K)) and the dis- 
pe rsion of the radial velocit y measurements {a{Vr)) given 
by iNordstrom et al.l l)2004j) to exclude as much as pos- 
sible bad quality radial velocity measurements or multi- 
ple systems out of Figure [S] In this sense, objects with 
e{Vr) > 1 km s~^ or with a{Vr) > 3 x e{Vr) were excluded. 
The remaining known binaries not excluded by these cri- 
teria are overploted by circles and triangles, respectively 
for FEROS and coude observations. Thus, excluding the 
known binaries and one point whose \Vrour — VfNordl > 
5 km s^^ we have a systematic difference 700 m s~^ and a 
r.m.s around the one-to-one relation of 800 m s~^ for the 
FEROS radial velocities. For the coude measurements the 
mean difference is about 1.4 km s^^ and the r.m.s about 
1 km s~^. These results show the very good agreement of 
ou r radial velocities with t hose obtained with CORAVEL 
bv lNordstrom et all l|20n4) and also that the FEROS and 
coude measurements can be combined since the systematic 
shift between the two instruments is within the errors. 

As a byproduct of the radial velocity measurements, 
projected rotational velocities V sm{i) were co mputed for 
the F EROS spectra using the calibration of iMelo et alJ 
l|200lh . For fast rotators {V sin i ^ 30 k m s~^) the fi - 
nal rotational velocities were derived as in iMeld l|2003|) . 
In Figure [5] (bottom) we compare our rotational veloc- 
ities with those from Nordstrom et al. (2004). Multiple 
stars are shown as encircled symbols. Again, the one-to- 
one relation is shown as a continuous line whereas the 
dashed-dotted lines mark the V sm{i) errors of 2 km s^^ 
for V sin j < 20 km s^^ and 10% of the observed ^ sin(i) 
for Vsin i > 20 km s~^. A closer investigation of the 6 
deviant points (those several km s^^ away from the one- 
to- one relation) show that based on the high a(K) given 
bv lNordstrom et all l)2004l) . these objects are Hkely to be 
multiple stars seen in blend. 



4.2. Spectral type, Li and Ha 

Spectral classification was carried out for the observed 
spectra using a combina tion of three d i fferent procedures: 
i) comparison with the Monies et al.l lll997 *) spectral li- 
brary, ii) calibrations of iTripic chio et al J1^97, 1999) for 
the Nal D and KI 769 9 A lines and, iii) the method devel- 
oped bv lTorrej 1 19991) in the Ha region. The good quality 
of the obtained spectral types can be seen in Figure 13 
where we plot the {V — I)c against our spectral types for 
dwarf stars. We estimate the error in our spectral type 
classifcation to be about one sub-class. Part of the spread 
in Figure |3| is real, as we are dealing with young stars, 
having distinct color excess, due to cool stellar spots, cir- 
cumstellar material, etc. 

Equivalent widths of the Li 6708 A and the Ha lines 
were measured from the FEROS or coude spectra using 
IRAF splot task. We estimate for low rotators {V sin i < 




Fig. 3. Our spectral types as a function of the {V — I)c- 



700 I I I I 



600 - 




500 - 



.400 - 



3 300 \- 



200 



100 



(V-I), 

Fig. 4. Distribution of Li equivalent width as a function 
of the {V — I)c- The lines represents the upper and lower 
limits for Pleiades members. 



20 km s ^) observed with FEROS an error of < 0.01 A and 
twice for stars observed with the coude. The distribution 
of the Li equivalent width as a function of the {V — I)c 
is given in Fig. 01 The large spread seen in Fig. 01 is ex- 
pected and it reflects the spread in age of the X-ray stars 
ffleuhau ser 19 97) . The combination of the Li equivalent 
widths and the color indexes was used to discriminate the 
sub-samples of young stars, the Pleiades age stars and the 
old stars, as can be seen in the figure. 

But in the search for young associations wc will con- 
sider "young star" if the Li equivalent width is at least 
90% of that observed for the Pleiades members for a given 
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{V — I)c- The few stars later than MO with undetectable 
Li were included in this young sub-sample as the Li may 
be co mpletely depleted fo r M stars already in the PMS 
phase l|Martin et al.lll994j) . 

In summary, from a sample of 1626 stars (1511 stars 
observed by us and 115 stars taken from the literature), we 
classified 255 as giants stars and 1371 as dwarf stars. We 
considered as dwarfs also the subgiant stars (many of the 
young stars are spectroscopically classified as subgiants), 
and 22 stars for which we did not obtain a luminosity 
class (mainly fast rotators or young stars taken from the 
literature). Among these 1371 dwarfs, 565, 282 and 524 
are considered respectively younger, having similar age or 
older than the Pleiades. 

5. Finding associations: the kinematical method 

Intuitively one can consider an association as a group 
of stars appearing concentrated together in a small area 
in the sky sharing some common properties such as age, 
chemical composition, distance and kinematics. However, 
if such a group is close enough to the Sun, its members will 
appear covering a large extent in the sky (as an example, 
Orion at 50 pc would cover almost the whole sky). Thus, to 
find a group, projected spatial concentration (i.e., in terms 
of right ascension and declination only) and proper mo- 
tions may not be enough. A better criterion is to look for 
objects sharing similar space motions (UVW) all around 
the sky (U positive towards the Galactic center, V posi- 
tive in the direction of Galactic rotation) . To obtain space 
velocities, proper motions {fJ-a, fJ-s), radial velocities and 
parallaxes are needed. Since most of the young stars ob- 
served in the SACY have no parallaxes, a convergence 
method to dynamically estimate their distances has been 
developed. This method is described below and applied to 
the BPA in the next section. 

5.1. The kinematical convergence nriethod 

As an initial step before applying the convergence method, 
we search for an initial point (UVW)o in the UVW-space 
around which a concentration of objects is observed indi- 
cating the existence of a possible association. This can be 
done basically in three different ways: 

i) considering only stars with known parallaxes; ii) con- 
sidering all stars with and without measured parallaxes; 
iii) using reliable information from the literature. 

i) In the first case, visual inspection of the distribution 
of objects in the UVW-space can already reveal cluster- 
ing around some points of the UVW-space. These points 
may furnish a starting point, (UVW)o, for the convergence 
method. 

ii) In the second case, all stars are considered regard- 
less whether they have parallaxes or not. For those with- 
out parallaxes a distance is found assuming a common 
age for all objects. (UVW) velocity components are then 
computed for all stars and as above, any observed concen- 
trations in UVW-space can be used as starting point. 



iii) A third source for the (UVW)o may be the liter- 
ature. For instance, the starting point for the BPA, pre- 
sented in Section 6, can be found either apply i ng (i) or 
taking the (UVW) given bv lZuckerman fc Sond l|2004l) . 

An additional ingredient needed to start the conver- 
gence method is the age which is estimated by the Li 
6708 A strength and the isochrone fitting. Once we have 
an (UVW)o and an age assigned we loop over all young 
stars of our sample. For each star we search the distance d 
which minimizes the quantity F{mv, Ha, fJ'S, Vr). The func- 
tion F is a quadratic combination of the evolutionary dis- 
tance (i.e., distance in magnitudes of a given star to the 
adopted isochrone) and the kinematic distance which is 
the Cartesian distance in km of the (UVW) of a star 
to the starting point (UVW)o: 



F(mt,,^a,^5,K.;d) 
[p X {M, 



" -I- 

(U ~ Uof + {V - Vof + {W - 



(1) 



Wo) 



211/2 



where p is a constant weighting the importance of the evo- 
lutionary distance with respect to the kinematic distance. 
The kinematic distance is, in other words, the modulus of 
the velocity vector 

dk = \i^-i^o I = [{U - Uo)' + (y~ Vo)^ + iW- Wo)']'/\2) 

A limit of 3 km s^^ (i.e., — fo| < 3 km s~^) corre- 
sponding to 3 times the typical velocity dispersion of open 
clusters is set as kinematical criterion along with F < 3.5 
which corresponds approximately to an evolutionary dis- 
tance of less than 2 mag. If both conditions are observed 
the star is classified as possible member. 

The procedure is then repeated until the convergence 
in the mean (UVW) is reached and the member list does 
not change any further. 

As the main concentrations in UVW-space are also 
spatially localized, albeit some of them covering large ar- 
eas in the sky, the distribution of the objects in physical 
space (XYZ) is used as an additional constraint. In some 
cases similar concentrations in (UVW) split up completely 
in the XYZ-space (Torre s et al. 2003). 

In general, ad-hoc limits in (XYZ) have to be imposed 
to avoid spurious stars and, for nearby associations in 
UVWXYZ-space, mutual contaminations. In these cases, 
the procedure is repeated with the new spatial constraints. 
A summary of the description given above is schematically 
presented as a flow-chart in Figure |5l 

It is worth noticing that due to the way our starting 
sample was defined (c.f.. Sec. 2) the members of the asso- 
ciations found by the SACY are a sub-sample (or a lower 
limit) of the full list of members. 

5.2. Membership probabilities 

We apply an alternative, objective, method to assess mem- 
bership probabilities of the stars in associations suggested 
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UVW taken 
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in the catalogue 




Compute new 
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Ad-hoc 
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Fig. 5. The convergence method. Left. As a first step, a concentration around a point in the UVW-space is found. 
The position of the stars defining this concentration at the HR diagram and their Li equivalent width give a first 
estimation for the age of the association. Right. Starting with an age and an initial (UVW)o determined in step 1, we 
look for the stars in the young sample which fit the association. 



by the kinematical convergence analysis. We employ a pre- 
dictive classification scheme that assumes group member- 
ship a priori for a certain sub-sample of objects, the so- 
called training set. For example, we can assume a training 
set derived by the kinematical convergence method. 

Our method works as follows: For each object, we stan- 
dardize all 6 variables UVWXYZ by subtracting a mean 
value and dividing by a standard deviation in order to 
guarantee similar weighting of each component in the 6- 
dimensional space. Means and standard deviations are de- 
rived from the training set only. Next, we apply a distri- 
bution free k-NN (k nearest neighbors) analysis to the 
entire sample to assign probabilities for each object. This 
approach dispenses with the need to introduce a probabil- 
ity density function, k — NN probabilities are calculated 
by defining euclidean spheres around each object in the 6 
dimensional space that contains exactly k objects. These 
k objects contain nt objects from the training set, and 
(trivially) k ~ nt other objects. The membership proba- 
bility a is then simply defined as a = nt/k. a is defined 
for each point in the 6-dimensional UVWXYZ-space, and 
changes discontinuously in steps of 1/fc. This method is, 
for example, used to analyze the large-scale spatial dis- 
tribution of X-ray s elected TTS from the ROSAT survey 
llSterzik et al.lll995l). and more details can be found in 
iMurtagh fc Heckl lll98/1^ . 

The use of the convergence method to define a bona- 
fide list of members is independent of the assessment of 
membership probability described here. Therefore, poten- 
tial members suggested by other works whose radial ve- 
locity and/or parallaxes have been collected from the lit- 



erature might be tested in similar way to those selected 
by the convergence method. 



6. Results: the (3 Pic Association 

As a test case for the convergence method described in 
the previous section, we look for the known BPA. As a 
starting point, only young stars with trigonometric paral- 
laxes are taken into account. Their UVW velocities show 
a concentration at U=-10 km s^^, V=-16 km s^^, W=- 
9 km s~^ which is very close to the space motion s found 
for the BPA defined bv lZuckerman fc Sond l|2004j) . 

A few proposed BPA members in IZuckerman fc Sond 
(2004) not a priori in the SACY sample (mainly A-F stars 
excluded predominantly due to the temperature criterion) 
were also included in our present analysis in order to test 
the consistence of the association as a whole, i.e., early- 
and late-type memb ers. In addition, membersh ip of a few 
ad-hoc stars not in IZuckerman fc SoTigl l)2004l) or in the 
SACY sample were also tested, for example, the case of 
the isolated cTTS V4046 Sgr. For all faint members of 
visual binaries, even wide, with no Hipparcos or Tycho-2 
data we used the proper motions and/or parallaxes of the 
primary star. 

In the convergence method we do not use the 
Hipparcos parallaxes when the errors are larger than 
2 mas, as this may indicate unreliable parallaxes. Anyway, 
these stars coud be used to check the resulting kinematical 
parallax. 
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6.1. Observational isochrone 



According to 'Zuckerm an et all ll200l£ 
12^4 Myr old, and Fcigcls on et'al] l|200e 



the BPA is 
estimate the age 
of the BPA member GJ 3305 to be 13^3 Myr. This is con- 
sistent with t he 11 Myr ob tained by the dynamical tracing 
back models (|Ortega et al. 2002. 2004). 

The col or-magnitude diag ram with the isochrone of 
10 Myr of ISiess et all l(2000l i is shown in Figure El In 
contrast to low mass stars which agree with the posi- 
tion of the isochrone, stars earlier than about G5 are 
up to 2 magnitudes off from the theor etical predictions. 
This di screpancy was already noticed bv lZuckerman et al.l 
^l2001a^ who overcame the problem by adopting observa- 
tional isochrones based on stellar populations with well 
defined ages. If on one hand isochrones are an important 
element to the convergence method applied here, enter- 
ing in the estimation of the evolutionary distance (see 
equation^, on the other hand, the temperature region 
where the differences between theory and observations 
seems more severe, is mainly beyond the temperature cut- 
off applied to the SACY sample. 

In spite of that, in order to be able to use the 
convergence method in future works covering a broader 
temperature ran g e, we adopted a similar strategy as 
IZuckerman et all l)2001a(l and use the BPA members of 
Table El to define an ad-hoc isochrone of around 15Myr. 




Fig. 7. Correlation between the X-position and U- 
component of space velocity vector showing an expansion 
of the members of the /3 Pic Association . Objects whose 
parallaxes are from Hipparcos Catalogue are overplotted 
with squares. 




(V-Ic) 



Our best fit gives (shown as a thick dot-dashed line in 
Figure EJ : 

M„ = 1.50+4.98(y-/)c-1.16(y-/)c+0.193(y-/)^(3) 

for -0.1 < {V -I)c < 3.1 

Eq. El is then the isochrone used in the convergence 
method. At a first glimpse, the way this observational 
isochrone was defined seems to lead to a circular argu- 
ment, namely, the use Eq. El into the convergence method 
to find the members which are going in turn to be used 
to refine Eq. El It is worth noticing that the objects used 
to set the early-type side of the observational isochrone 
all have determined parallaxes and are anchored in the 
UVW space and on the color-magnitude diagram. The 
intrinsic assumption is that these early-type objects are 
indeed members of BPA. According to our membership 
probabiliti es given in Table El all bu t one A and F stars 
taken from lZuckerman fc Sond l)2004(l have a membership 
probabilities higher than 85% and are considered to be 
bona-fide members which validates our approach to define 
the observational isochrone. 



Fig. 6. HR diagram from the proposed members of the 
P Pi c Association. Th e isochrones over-plotted are those 
from lSiess et al ]f2000). The dotted and dashed lines cor- 
respond to the isochrones for 5 Myr and 50 Myr, respec- 
tively. The isochrone for 10 Myr is shown as thick solid 
line whereas the observational isochrone is represented as 
a thick dotted-dashed line. 



6.2. Evidence for expansion 

A plot of the velocity component U of the space velocity 
vector against the space direction X is seen in Figure d 
The farthest components show a larger U-component ve- 
locity suggesting that the association is actually expand- 
ing. A quick check of whether this effect is an artifact of 
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Fig. 8. Combinations of the sub-spaces of the UVWXYZ-space showing a well defined clustering in both kinematical 
and spatial coordinates. 



our convergence method can be carried out looking for the 
objects whose parallaxes were taken from Hipparcos and 
not derived by the convergence method. These objects are 
overplotted with squares in Figure It is clear that the 
correlation does exist, even if the members are restricted 
to 50 pc as is the case of the IZuckerman fc Son3 l)2004|) 
analysis. A linear least-squares fit gives: 

C/ = 0.053(A:) - II.l (4) 

Such a correlation seems to be a common trend ob- 
served among the young stars (i.e., age ^ agepieiades) of 
the SACY sample in contrast to the older sample where 
no correlation is seen. In addition, no significant correla- 
tion seems to exist either between Y and V or W and Z, 
regardless of the age of the objects. Interpretation of what 
caused the trend seen in Figure [3 and its possible links to 
star formation and early dynamic interaction is postponed 
for a future paper. 

In the convergence method described in Sec. 5.1, Eq.^ 
had to be taken into account regardless the physical origin 
of this expansion. This is done simply by replacing the Uq 
in Eq.Hlby an Uq such that UQ = Uo + kX. 



6.3. The (3 Pic Association membership probabilities 

We estimated the membership probabilities of all possible 
members to the BPA as discussed in Sec. 5.2. We assume 
the training set as derived by the kinematical convergence 
method, given in Table 

For the BPA analysis we assume k = 20. All mem- 
bers of the training set have probabilities in the range of 
0.7-1.0, i.e. they are already well confined and clustered 
in the UVWXYZ-space. TYms^hona-fide members in this 
list have membership probability of, at least, 0.75. These 
probabilities are displayed in the last column in Table 

We do find four stars above this probability thresh- 
old which are not given in Table Three of them have 
magnitude differences with respect to the observational 
isochrone (Eq. |2l larger than 1.4 mag and therefore were 
discarded. The other is Hip 88726, having (X, U) values 
(42.3pc, -13.0 km s~^), not compatible with the observed 
trend seen in Fig. (or Eq. 0J). We rejected it as a BPA 
member in spite of a mem bership probabilit y of 90% and 
an acceptable magnitude. ISong et all l)2003j) also consid- 
ered it as an outlier. 

The final list of high probability proposed members is 
the same as that obtained by the convergence method. 
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6.4. The SACY f3 Pic Association 

The final results for the BPA are summarized in Table El 
where the magnitudes and colors are corrected by the pres- 
ence of companions as explained in Section 2.2. 

Table contains 41 highly probable members, that 
could be used to define the main properties of the BPA 
an d 2 possible member s. 

IZuckerman fc Sonel ( 20041) proposed a list of 33 mem- 
bers'^ to the BPA including stars hotter than the SACY 
limit of B — V > 0.6 and in both hemispheres. Only one 
very faint star proposed as member by Zuckerman & Song 
was not used in our analysis as it lacks the kinematical 
data needed to be injected in the convergence method. 



^ Actually their list has 34 entries, but 2 entries are for both 
members of the SB2 V824 Ara. 



Table 6. Data of the proposed members of the /3 Pic Association. Column 1: Star Identification; Column 2 and 3: Equatorial Coordinates; Column 4 and 5: 
Proper Motions; Column 6: Radial Velocities and Notes about duplicity; Column 7: Rotational Velocities; Column 8 and 9: Photometry; Column 10: Spectral 
Types; Column 11: Li Equivalent Widths; Column 12: Trigonometric Parallaxes; Column 13: Kinematical Parallaxes; Column 14 to 16: Space Velocities; Column 
17 to 19: Space Positions; Column 20: Function F (see text); Column 21: Membership Probabilities; Column 22: Notes. 
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Hemis phere , and not observe d by us ; S = observed in the SACY p r ogram ; H = photomet r y from Hip parcos; N = rad i al velo city and /o r V sin i from [Nords trom et al.1 
ExMh: rn IZuckerman et a.l .1 (l2nm ah : ^2^ iRa.rbier-Rrossa.t Fi^ tmh: (^) ISon^et a.l l E(K)^: (4) icutisnotn et a.lJ Hmh: IWei^ (Il991li: ffi^ lla.nn j il 97qh: (7\ 
[ Robertson fc HamihonI lll987h: fSllKoen et al.l J2002D (9) IWeij Jl991ah: (10) ICizis et alJ mm: (111 iFeigelson et aO teOOeT): fl2) IChristian fc Mathioudakig tOQ^: (13 ) 
iFavata et al.' figgS*); (14) 'Wichmann et al.' (200f): (15) Strassmeier & Rice ('200 d|); (16)lQuast et all ('2000'): fl7l 'Erspamer fc North' ('2003'): flSl lGrav fc GarrisonI ilQS'Al : 
ri9 :.CutisD0tQ (,1998.) : (20) (,Eggen...l97&) : ('21) .Cutispoto et al. (,1999 ): (22) Delfosse et al.l I^1999^ : dg '.CutispotQ. 1,1998) : (20) (,.Eggen.l97S. ). 



C. A. O. Torres et al.: SACY: Sample & Searching Method 



13 



Four stars of their list arc not proposed as probable 
members: 

- Hip 88726 and Hip 79881 - they do not fit into 
the expansion relation (Hip 79881 has a membership 
probability of 65%) and they were also rejected by 
ISong et all lE)031 

- HD 203 and HD 199143 - they arc high rotational ve- 
locity F stars (T^sin(i) = 155 km s~^) and their radial 
velocities are less reliable. Both stars would be con- 
sidered as good members with radial velocity values 
within one sigma of the mean used values or if we have 
used published radial velocities. We consider them as 
possible members. 

From the 13 new members proposed in this paper only 
two have Hipparcos parallaxes (HD 168210 and Hip 23200 
= V1005 Ori). Another one is the new visual companion 
of HD 164249, and we use the Hipparcos parallax of the 
primary. 

In Figure |H1 we show different projections in 
UVWXYZ-space. In all panels a well defined clustering 
in both spatial and kinematical coordinates is seen which 
gives confidence of the solution obtained by the conver- 
gence method. The clustering has an oblate shape along 
the X and U directions, as discussed in section 6.2. 

According to Table El the association has a mean 
(UVW) of (-10.1 ± 2.1 km s-\ -15.9 ± 0.8 km s^^, 
-9.1 ± 0.9 km s-i), and a mean (XYZ) of (20 ± 34 pc, 
—5 ± 15 pc, — 15 ± 7 pc). It spreads over a 100 pc in- 
terval in X, 50 pc in Y and 30 pc in Z. The distances 
of the proposed members of the BPA to the Sun range 
from 10 pc up to 80 pc. If only stars later than GO 
(28 stars) are considered the mean (UVW) would barely 
change (-10.1 ± 2.0 km s-\ -16.0 ± 0.7 km s'^ and 
—9.1 ± 1.0 km s^^). The distances would change less than 
1%. 

It is evident that the solar position is almost inside 
the BPA and thus it could have members spread over all 
the sky, and there may exist members undetected by our 
survey due to its partial sky coverage. Nevertheless, the 
distribution of the members in XZ-plane seems to indicate 
that the BPA is an association mainly in the Southern 
skies. Even a very high negative X value, if we extrapo- 
late the expansion using the mean Y and Z, would imply 
in stars near a ^ 05H and 5 ^ + 10°. To be sure 
where BPA ends at negative X direction it is essential to 
extend SACY at this region. The gap near X = -10 in 
Figure comes from the incomplete sky coverage in this 
region. Nevertheless, an extreme positive X value would 
fall within our sample. 

Figure El shows the Li equivalent width as a function 
of color index for all proposed members of Table O As 
seen in Figure El they lie above the upper envelope set by 
the Pleiades members, as expected from the fact that by 
definition the young sample contains only stars showing a 
considerable Li. Nevertheless, it is worth noticing that for 
the G-K dwarfs (where the Li-age correlation is stronger) 
the observed equivalent widths are much higher than the 
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Fig. 9. Distribution of Li equivalent width as a function 
of the {V — I)c for the proposed members of the /3 Pic 
Association shown in Table El The hatched area shows 
the Li equivalent width range observed for the Pleiades 
members. 

Pleiades upper envelope indicating a much younger age for 
these objects relative to the Pleiades. We tested the con- 
vergence method including the Pleiades age sample with- 
out finding new members. 

To obtain the Li abundances, presented in Table [71 we 
used our modified version of the known LTE code made 
by Monique Spite, of Paris-Meudon Observatory. These 
Li abundances were determined by comparing the equiv- 
alent widths of the synthetic lines to the measured ones, 
on the observed spectra. The method is very simple: given 
an equivalent width for a given line, we calculate its the- 
oretical profile for a given element abundance and its cor- 
responding equivalent width. If the theoretical equivalent 
width is not good, the process is retrieved for an interpo- 
lated abundance value, until the difference between both 
values be acceptable. Our criterion is that this difference 
be less than 0.2 mA. In the computation of the synthetic 
profile we take in account the 4 components of the Li 
resonance line, bei ng the wavelengths ari d the oscillators 
strengths given by ijAndersen et al.l ll984V The used stel- 
lar atm ospheric models are from Gustafsson and collab- 
orators ijEdvardsson et al1ll993(l . The effective tempera- 
tures used for the atmosphere models were determined 
from {V — I)c color index using the (Kenvon Sz HartmannI 
I199,'tI) calibration. The other atmosphere model parame- 
ters were fixed a priori: the metallicity as [Fe/H] = 0.1 
( Castilho et aL, 2005) and log g = 4.5, for stars of lumi- 
nosity class V, or log g — 4.0 for luminosity class IV. 
The microturbulence velocity also was fixed as 2 km s~^. 
The Li abundance is not very sensitive to gravity and the 
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Fig. 10. Distribution of Li abundances as a function of 
the effective temperature for the proposed members of the 
f] Pic Association (Table[7|). The hatched area shows the Li 
abundances with range observed for the Pleiades members 

main source of errors is the inaccuracy of the equivalent 
widths (specially for fast rotating stars) and the effective 
temperatures. 

A further evidence for the youth of the BPA members 
is given by the projected rotational velocities {Vsm{i)). 
They span over a large range from ^ 5 up to 150 km s~^ 
as typically observed in yo ung clusters of similar age such 
as IC 2391 and IC 2602 llStauffer et al.lll997tl and the 
slightly older a Per (lBouvier^ll997^ . This spread is also in 
good agreement with the theoretical angular momentum 
evolution tracks and are the result of different time-scales 
for disk-unlocking which is probably related to the disk 
life-time itself (e.g. ..Allain et al.■.1997^) . 

6.5. Comments about some proposed members 

AT Mic and AU M ic were already classified as young 
stars by EggerJ l)l968|) whereas their connection w ith Pic 
was first proposed by iBarrado v Navasciies et al. (J^SS) . 
Errors in the Hipparcos distances for the AT Mic pair are 
very high and their proper motions have differences larger 
than usual between Hipparcos, Tycho-2 and UCAC2. The 
pair has an orbital solution with a period of 209 yr in the 
Sixth Catalog of Orbits of Visual Binary Stars (Hartkopf 
& Mason, htpp://ad. usno.navy.mil/wds/orb6. html), thus 
catalogs with larger time base may have the proper mo- 
tions contaminated by the orbital motion. Therefore, we 
prefer to use the proper motions of Hipparcos, as they bet- 
ter reflect the center of mass motion. Our kinematical dis- 
tance is close to that of the companion AU Mic. The weak 
Li line of AU Mic was first suspected bv lde La Reza et all 



Table 7. Li abundances, in the scale log A(H)=12, for 
the BPA proposed members having EW^i > in Table E| 
Column 1: Star identification; Column 2: Effective tem- 
perature, from {V — I)c', Column 3: Equivalent width of 
the resonance lithium line in mA; Column 4: log of Li 
abundance. 
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in an attempt to explain the odd Li-rich red dwarf 
V1005 Ori (Hip 23200; Gl 182) by means of Li produc- 
tion by spallation reactions. They discarded this possibil- 
ity and proposed that V1005 Ori may be a member of a 
very young kinematical group. These stars may now have 
their oddities explained as proposed members of the young 
BPA. 

WW PsA and TX PsA form a low mass wide visual 
binary, the first one being slightly brighter and hotter. 
TX PsA has a high Li equivalent width (the point around 
(U - I)c 3 in FigurcEl. In contrast, WW Psa has the 
Li fully depleted. Song et al. (2002) interpreted the pair in 
the context of Li depletion boundary and concluded that it 
can give a strong observational constraint to age determi- 
nations using the Li depletion boundary. Our kinematical 
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distances put the pair closer to the Sun than the Hipparcos 
parallax which has a large uncertainty. If our distance is 
correct, the stars could be fainter, and this should be con- 
sidered in the discussion of the Li age determinations. The 
faint companion of HD 139084 is very similar to TX PsA 
(the spectral type may be somewhat lat er) and it can help 
in this context. The {V - I)c given bv lSong et all l)2003(l 
seems too blue; from 2Mass, the two stars have similar 
(V-K), and we use in the paper the TX PsA color for it. 

HD 199143 (Hip 10331 1) and AZ C ap (HD 358623) 
were recently proposed by iKaisler et al.l (|2004) as mem- 
bers of the BPA. HD 199143 is a fast rotator and its radial 
ve locity has a large spread. If w e had used the value given 
bv lvan den Ancker et al.1 l)2000tl . the results yielded by the 
convergence method would support the BPA membership 
for both HD 199143 and AZ Cap. Our results challenge 
the existence of the propose d Capricornius moving group 
ijvan den Ancker et al.ll200o|) . which is more likely a sub- 
sample of the BPA. 

•q Tel is a very fast rotator. The V smji) in Table El was 
estimated by ICrav & CarrisonI ( 1987f) . Thus its radial ve- 
locity is unreliabl e and we used km s~^ as proposed by 
ISong et al]l)2003|) . This value is very close to the radial ve- 
locity of its wide companion HD 18132 7. It has a propos ed 
brown dwarf companion at 4.1" (Lowrance et ahlboOOl) . 

There are 5 double line spectroscopic binaries (SB2) in 
the BPA: 

Hip 23418 was found bv lDelfosse et all l)l999|) to be a 
SB2 with a period near 12 days and mass ratio of 0.57. 
In Table we used their systemic velocity. The star is 
also a very close visual double detected by Hipparcos, the 
secondary star being one magnitude fainter. In TableElwe 
try to give the magnitude and color of the primary for this 
triple system. 

HD 35850 was found by iNordstrom et all l(2004l) to 
be a SB2, who give a mass ratio of 0.72, the systemic 
velocity, but no rot ational velocity. For V sm.(i) we use 
that determined bv IWich mann et al ] ll2003t) . The star is 
considered an outlier bv lSong eT^ J(|2003^ since they did 
not know about the radial velocity variati ons. It is also a 
variable star with a peri od near one day l|Budding et al.1 
I2OO3': 'Wright et al."2004h. 

HD 155555 (V824 Ara) is actually a multiple sys- 
tem that includes a well studied SB2 (G7+K0IV) 
fst rassmeier fc Ricell200nl) an d a fainter optical compan - 
ion, proposed as member bv IZuckerman fc Sond ((20041) . 
We observed also this fainter companion, a M4Ve star 
with no detectable Li line, and our radial velocity iscon- 
sistent with the value obtained bv IPasauini et al ] lll99lh 
(3.5 km s~^). In Table we used for its proper motions 
and parallax the data of the brighter component. 

V4046 Sgr is not in the SACY sample, since it is not in 
IRXS, but we included it to test if it could belong to some 
association as it is an isolated classical TTS. It is also a 
SB2 with a circu mbinary disk (IStemoels fc Gahml l2004i 
)uast et al.ll200o() . The distance obtained bv lOuast et alJ 
l|200nl) . 83 ± 8 pc, evolutionary model dependent, agrees 
well with our kinematical one of 72 pc. The age found by 



IStempels fc Gahml l|2004|) also agrees well with the age of 
the BPA. GSC 7396-0759, a nearby young red star and op- 
tical counterpart of the ROSAT source IRXS J181422.6- 
324622, could be a distant (2.8') physical companion of the 
V4046 Sgr system, based on its radial velocity, Li equiv- 
alent width and magnitudes. Using for it the V4046 Sgr 
Tycho-2 proper motions, the star agrees well with a BPA 
membership. However, the shape of its CCF suggests it 
may be a SB2 and, in this case, it could not be a BPA 
member. 

HD 161460 is a SB2 (KOIV + KlIV), observed only 
twice, one near conjunction and in the other spectrum 
the lines are not very well resolved, but both Li lines are 
strong. The radial velocity is reliable and the mass relation 
is near 0.5, implying a later spectral type for the secondary 
star. 

BD+ 05 378 is a single line spectroscopic binary de- 
tected bv lSong et al.l l|2003() whose estimated systemic ra- 
dial velocity is used here. We observed it only once, after 
its proposed inclusion in the BPA, and the other spectro- 
scopic information in Table E| comes from our spectrum. 

Thus, there are 6 known spectroscopic binaries, 6 close 
visual binaries (separation less tha n 200 AU) and 10 
wide v isual binaries (see Table and IZuckerman fc Sond 
l|2004) ) in BPA over 41 (or 43) proposed members. 

7. Summary and conclusions 

In the present paper we present the SACY, a high- 
resolution optical spectroscopic survey aiming to look for 
undiscovered nearby associations among the optical coun- 
terparts of the ROSAT X-ray sources in the Southern 
Hemisphere. From an initial list containing 9574 IRXS 
sources in the Southern Hemisphere, we find 1953 of them 
having counterparts with B — V > 0.6 in the Tycho-2 or 
Hipparcos Catalogues. During the last seven years, 1511 
objects were observed at least once. Published information 
for 115 objects was used to complete our measurements. 
SACY is ~ 80% complete in the Southern Hemisphere. 
The main results of this paper are summarized below. 

— Radial velocities and Vsai(i). Accurate radial veloci- 
ties {(j{Vr) ^ 0.3 km s~^) have been derived for all 
observed sample using FEROS and/or coude observa- 
tions. Good quality V sm{i) were derived from FEROS 
observations. Both quantities, radial velocities and 
ysin(i), were compared to recent measuremen ts ob- 
tained with CORAVEL ((Nordstrom et alJl2004^ show- 
ing a very good agreement. 

— Spectral classification. Spectral type classification was 
carried out using a combinatio n of thr ee different 
methods: i) comparison with the iMontes et al. ( 199^ 
spectral lib rary, ii) calibrations of I'Mpicchio et alJ 
(Il997tll999l) for the Nal D a nd K I 7699 A hnes and, 
iii) the method developed by iTorresI (|1999.) in the Ha 
region. 

— Young stars. Based on our spectroscopic observations 
Li equivalent width could be measured for the first 
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time for most of the 1256 dwarfs observed in the survey. 
Using also the information for 115 stars from the liter- 
ature and based on the Li criterion we consider, 565, 
282 and 524 dwarf stars as respectively younger, hav- 
ing a similar age to or being older than the Pleiades. 

Since parallaxes exist only for a small part of the sam- 
ple, we developed a kinematical convergence method to 
overcome this problem allowing us to derive dynamic par- 
allaxes. As an example, this method is applied to an initial 
concentration of young stars at U,V,W=-10,-16,-9 km s"^ 
which i s very close to the propos ed mean velocities for 
BPA bv lZuckerman fc Sond l)2004() . 

— Defining the BPA in (UVW). According to TableEl the 
association has a mean (UVW) of (— 10.1±2.1 km s~^, 
-15.9 ± 0.8 km s-\ -9.1 ± 0.9 km s-^). 

— New members of the BPA. Our defin i tion o f the BPA is 
close to that of IZuckerman fc Sond (|20^. Using our 
convergence method and data from the literature we 
verify that most of their northern and early-type mem- 
bers fit well into our definition of the BPA. According 
to our results we propose 13 additional members, in- 
cluding the well known active stars V1005 Ori and 
V4046 Sgr. This latter is not a "SACY star" and was 
formally considered as typical isolated TTS. Our re- 
sults suggest that V4046 Sgr belongs indeed to the 
BPA, being one of the oldest cTTS known up to now. 
We propose a wide companion red dwarf star for it. 
We find also a faint companion for HD 164249. We 
propose a list of 41 bona-fide BPA members. 

— Observation al isochrone of th e BPA. The 10 Myr 
isochrone of ISiess et all l|200nt) fits part of the tem- 
perature range of the members of the BPA. We then 
reversed the process and used the color-magnitude di- 
agram of the BPA t o define an ad-hoc obser vational 
isochrone as done bv lZuckerman et alJ l)2001a|) . 

— Indication for expansion. We find a positive correla- 
tion between the X-direction and the U-component of 
the space velocity vector for young stars in the SACY 
sample. This behavior is most clearly seen in the case 
of the BPA. Analysis of this expansion and its impli- 
cations will be discussed in a future paper. 

— Membership probability. We have developed a mem- 
bership assessment procedure based on a predictive 
classification scheme to be used with the convergence 
method. Both the value of the F and membership 
probability agree well. 
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